Self-generated magnetic flux in YBa 2 Cu 3 07_ x grain boundaries 
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Grain boundaries in YBa2Cu307-i superconducting films are considered as Josephson junctions 
with a critical current density j c (x) alternating along the junction. A self-generated magnetic flux is 
treated both analytically and numerically for an almost periodic distribution of j c (x). We obtained 
a magnetic flux-pattern similar to the one which was recently observed experimentally. 



I. INTRODUCTION 

Grain boundaries in high-T c cuprates are interesting 
and important for both fundamental physics and applica- 
tions of high-temperature superconductivity.B~H Conven- 
tional models of strongly coupled Josephson junctions are 
applicable to describe electromagnetic properties of fflffin 
boundaries in thin films of high-T c superconductors OU A 
remarkable exception of this rule is the [001]-tilt bound- 
ary in YBaaCjUaOy-^ films with a misorientation angle 
close to 45 .ull3 Indeed, these grain boundaries have an 
anomalous dependence of the critical current I c on an ap- 
plied magnetic field H a . Contrary to a usual Fraunhofcr- 
type dependence I c (H a ) with a major central peak at 
H a = and minor symmetric side-peaks the asymmet- 
ric 45° [001]-tilt grain boundaries demonstrate a pattern 
without a central dominant peak. Instead two symmet- 
ric major side-peaks appjear at certain applied magnetic 
fields H a = ± H sp ^ oWb 

Several mechanisms have been suggested to explain 
this phenomena.ffiM The anomalous dependence I c {H a ) 
with symmetric major side-peaks is obviously a result 
of a specific heterogeneity of electrical properties of the 
asymmetric 45° [001]-tilt grain boundaries. Two funda- 
mental experimental observations in..conjunction explain 
this heterogeneity in a natural wayJjll3 First, a fine scale 
faceting of grain boundaries was discovered in experi- 
ments usina-thj&iransmission electron microscopy (TEM) 
technique BOlIa These facets have a typical length-scale 
I of the order of 10-100 nm and a wide variety of orien- 
tations relative to the axis of symmetry of the supercon- 
ductor. Second, quite a few of recent experiments provide 
an evidence of a predominant d x 2_ y 2 wave symmetry of 
the order parameter in many of the high-T c cuprates. In 
some experimental studies the symmetry of the order pa- 
rameter is more complicated and is shown to be-a-cectain 



mixture of the d x 2_ y 2 and s wave components. cE3~c3 

These two fundamental experimental observations in- 
dicate the existence of two contributions to the phase dif- 
ference of the order parameter across the grain boundary. 
Indeed, consider a meandering grain boundary in a film of 
a superconductor with the d x 2_ y 2 wave symmetry of the 
order parameter and assume that there is a certain mag- 
netic flux inside the grain boundary. In this case there is a 
phase difference ip caused by the magnetic flux and at the 
same time there is an additional phase difference a caused 



by the misalignment of the anisotropic d x 2_ y 2 wave su- 
perconductors. Therefore, the tunneling current density 
j c is defined by the total phase difference A = ip — a. 
A model describing this Josephson current density j(x) 
results from an assumption that j(x) oc sin [<p(x) — a(x)], 
where x is along the grain boundary lincJij The local val- 
ues of the phase difference a(x) depend on the relative 
orientation of neighboring facets. In the case of an asym- 
metric 45° grain boundary we hasse a(x) — or it, and 
therefore j(x) oc sin ip(x) cos a(a;).Ej In other words in the 
framework of a model relating j(x) to the orientation of 
the facets we arrive to j(x) — j c {x) sanp{x) with an al- 
ternating critical current density j c (x) oc cosa(a;). The 
dependence j c (%) is imposed by the sequence of facets 
along a grain boundary line. If this sequence is periodic 
or almost periodic then the function j c (x) is a periodic or 
almost periodic alternating function. The typical length- 
scale for j c (x) is of the same order as the length of the 
facets I, i.e., this length-scale is of 10-100 nm. 

Variation of orientation of facets along a meandering 
grain boundary leads to formation of local superconduct- 
ing current loops even in the absence of an applied mag- 
netic field if the total phase difference A ^ 0. It was 
predicted, in particular, that these current loops can gen- 
erate a certain magnetic flux at a contact of two facets 
with a = and a = tt.u 

Self-generated randomly distributed magnetic flux was 
discovered in asymmetric 45° [001]-tilt grain boundaries 
in YBa2Cu307_a; superconducting films in the absence 
of an applied external magnetic fieldH3 This flux cf) s (x) 
changes its sign randomly and has an amplitude of vari- 
ations less than the flux quantum c/jq. The average value 
of 4> s {x) along the grain boundary is nearly zero. Notice- 
able, that this disordered non-quantized magnetic flux 
was observed only for the samples exhibiting the anoma- 
lous dependence of the critical current I c (H a ) on mag- 
netic field with the two symmetric major side-peaks. 

It was shown analytically that under certain condi- 
tions a stationary state with a self-generated flux ex- 
ists for a Josephson junction with a—periodically alter- 
nating critical current density j c (x)E5 The same spatial 
distributions of j c (x) result in an anomalous dependence 
of the. critical current I c (H a ) on the applied magnetic 
fieldO Numerical calculations show that two symmet- 
ric major side-peaks appear for a periodically alternat- 
ing j c {x). The randomness of the critical current den- 
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sity j c {x) smears these peaks but leaves their position in 
place at weak randomness. We therefore conclude that 
the experimental observation of the -well pronounced ma- 
jor side-peaks on the curve I c (H a )ul!3 means that the 
alternating critical current density is a periodic or al- 
most periodic function of a;. A noticeable randoinness of 
j c (x) would smear out the dominant side-peaks. 

In this paper we calculate both analytically and numer- 
ically the self-generated flux <f> s in a Josephson junction 
with an almost periodically alternating critical current 
density j c (x). The paper is organized as follows. First, 
we review briefly— the case when j c (x) is a periodic alter- 
nating function.L3 We derive the main equations of the 
two-scale perturbation theory and apply these equations 
to analyze the non-quantized self-generated flux. This 
approach forms a basis for the following analytical calcu- 
lations. Next, we treat the self-generated flux for the case 
of an almost periodic alternating critical density j c (x) 
and start with a qualitative approach to the problem. 
We review then the results of our numerical simulations 
which verify the qualitative consideration and exhibit a 
magnetic flux-pattern which is similar to the one that 
was recently observed experimentally.! 2 ^ 

II. MAIN EQUATIONS 

It is convenient for the following analyses to write the 
function j c (x) in the form 

jo = {jc)[i + g(x)}, (l) 

where (j c ) is the average value of the critical current den- 
sity j c {x) over an interval with a length L 3> I 

(jc) = ^J o Jc(x)dx. (2) 

The function g(x) introduced in Eq. (|l]) alternates on a 
typical length scale of I. Note that by definition the aver- 
age value of g(x) is zero, i.e., (g(x)) = 0. The maximum 
value of \g(x)\ varies from \g(x)\ max > 1 to \g(x)\ max > 1. 
We assume also that A < / C Aj, where A is the London 
penetration depth and 

is an effective Josephson penetration depth. It is worth to 
mention that in the case of an alternating current density 
this effective penetration depth is not a local characteris- 
tics of a tunnel junction. It is rather an effective quantity 
defined on the same typical length scale as (j c ). 
The phase difference ip(x) satisfies the equation 

A 2 r ip" - [1 + g{x)] sin tp = 0. (4) 

In the limiting case I <C Aj it is convenient to write a so- 
lution of this equation as a sum of a certain smooth func- 
tion il>(x) with a length scale of order Aj and a rapidly 
oscillating function with a length scale of order Zcil 



(p(x)=il;(x)+t(x). (5) 

We assume also that \£(x)\ <C IVK 2 -)!- After substituting 
Eq. (||) into Eq. (|]) and keeping the terms up to the first 
order in £(x) we obtain 

AjiP" + A 2 j£" - [1 + g(x)] [sin ip + £ cos ip] = 0. (6) 

Note, that experimentally! 2 ! the self-generated flux was 
observed by a SQUID pickup loop with a size of several 
Aj ^> /. It means that this method is averaging out 
the fast alternating flux defined by the phase £(x) and 
measures the spatially smooth flux defined by ip(x). 

Next we consider briefly the case of a periodically al- 
ternating critical current density j c {x) which forms the 
basis of the following analysis of a general case with j c (x) 
being a randomly alternating function. 



III. PERIODICALLY ALTERNATING CRITICAL 
CURRENT DENSITY 

A. Two-scale perturbation theory 

If the critical current density j c (x) is a periodic func- 
tion, then g(x) also is a periodic function. In this case 
an approximate solution of Eq. (||) caa be obtained based 
on a two-scale perturbation theoryES As a first step in 
order to apply this approach to Eq. (||) we separate the 
fast alternating terms with a typical length scale I and 
the smooth terms varying with a typical length scale Aj 

( A 2 jiP" - sin ^ - g£ cos ip) + (A 2 jC' - g sin iji) = 0. (7) 

In Eq. (^) we omitted two out of three fast alternating 
terms of Eq. (JsJ) since they are proportional to £(x) and 
therefore are smaller than the term proportional to g{x). 
Next, we note that the terms included into the first pair 
of brackets in Eq. (fjj) cancel each other independently 
on the terms included into the second pair of brackets 
in Eq. (^) as these two type of terms have very different 
length scales / and Aj and I <C Aj.E3 The same reasoning 
is applicable to the terms included into the second pair of 
brackets in Eq. (f7j). As a result Jve obtain the following 
two equations for and ip(x)c3 

A 2 J £"=.g(x)sinV, (8) 
A 2 7 ?/>" - sin if) - (g(x)£(x)) cos ip = 0. (9) 

It is worth to note that we obtain the two functions ip(x) 
and from one equation ([?]) as only two type of terms 
with different typical length scales / and Aj appear in 
Eq. (|^). If g(x) would have a wide range of typical length 
scales the above separation would not be possible. 
Introducing the Fourier transform of g(x) as 
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we find the solution of Eq. (|J) in the form 



simp 
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-£ g (x) sinip, (11) 



where the sums in Eqs. ( |10P and ( |11[ ) are over the wave 
vectors k = 2im/£, C is the length of the junction, and 
n is an integer. It is worth mentioning that the func- 
tion £ g (x) is defined only by the alternating components 
of the critical current density j c (x). Also while deriving 
Eq. dll] ) we ignored the spatial dependence of sin?/'- This 
can be done since on the length-scale I the variation of 
the smooth function sinip(x) is of order l/Aj <C 1. The 
alternating part of the critical current density has the 
typical wave numbers k ~ l/l. Therefore, using Eq. (|ll| ) 
we estimate £(x) as 



£(x) smipj^gix) 



(12) 



It follows from this estimate that the typical values of 
the phase difference £(x) are small ((|£(x)|) <C 1) if 



(Iff(s)IX^. 



(13) 



Next, using Eq. ( ) we rewrite Eq. (g) for the smooth 
phase shift ip(x) in the form 



Aji/j" — sin ip + 7 sin ip cos ip = 0, 



where 



7 = (g(x)tg(x)) 



simp 



(14) 



(15) 



is a constant. A similar derivation for the current density 
across the tunnel junction results in 



j(x) = (j c ) sin -0(1 - 7 cos?/)). 



(16) 



The magnetic field B s (x) generated by the alternating 
component of the critical current (j c ) g(x) is given by 



(17) 



and the averaged field (B s (x)) = 0. The alternating mag- 
netic flux 4> s produced by the field B s is equal to 
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Combining Eqs. (|l5|) and ( |l7| ) we find for 7 the formula 



cX 



(B 2 S ) = 



(Bl) 



(19) 



where we introduce a characteristic magnetic field 



47T . 

Bj = — {jc Aj. 

c 



(20) 



It follows from equation (|l9| ) that 7 is a positive constant 
which can be estimated as 



(21) 



The energy of a Josephson junction £ takes the formE£l 
£ = £0 + wnere £0 is independent on tp(x) and 



(To) 
2e 



j dx[ l - A>' 2 - [1 + g(x)} cos ^} . (22) 



Using Eqs. (g), ([il|), and the definition of 7, we obtain 
the energy £ v in terms of the smooth phase shift ip(x) 



f — 



4e 



dx^A 2 jip' 2 - 2cos-0 -7 sin 2 ?/>}. (23) 



Note, that solutions 0(x) of Eq. ( |i4[ ) correspond to the 
minima and to the maxima of the energy functional (p3) . 



B. Non-quantized self-generated flux 

Let us apply Eqs. (g), (||), and dl5| ) to a consider the 
stationary states of a Josephson junction with a certain 
length C ^> Aj in an absence of applied magnetic field. 
In this case the average flux inside the junction is zero 
and thus an alternating self-generated flux S (x) appears 
simultaneously with a certain phase ip = const as it fol- 
lows from Eqs. (0) and @. 

In the stationary state with ip — const the values of ip 
are determined by Eq. (ph which takes the form 



simp (1 — jcosi/j) = 0. 



(24) 



Note, that this equation means also that the current den- 
sity j(x) across the tunnel junction is equal to zero. 

In the case 7 < 1 equation has two solutions, 
namely, ip = and ip = 7r and thus, as follows from 
Eq. (|il"|), there is no self-generated flux. It is also worth 
mentioning that the energy of a Josephson junction £ has 
a minimum for ip = and maximum for ip = 7T. 

In the case 7 > 1 there are four solutions of equation 
namely, ip = —tpy, 0, ipy, n, where 



ip 1 = arccos (1/7)- 



(25) 



The energy £[ip{x)] has a minimum for ip = iip 7 and a 
maximum for ip — 0, 7r. The self-generated flux 
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— ^smV> 7 = T<Po^r- (26) 
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arises in the two stationary states with tp = if/S-i each 
of these states corresponds to a minimum energy £. The 
assumption <C 1 restricts the value of 7. However, 

it follows from Eqs. © and (§|) that (£(a;)) < 1 and 
7 > 1 hold simultaneously only if 



A 7 ,. f \ 1 1 At 

-f < <ls(*)l> « -f- 



(27) 



Using equations ( |12| ) and ( |26| ) we estimate |</> s (a;)| as 

7 Z 2 



|0 s (a;)| ~ <p 
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2 < 0o- 



(28) 



The above results hold only for a periodic critical cur- 
rent density j c (x) and the predicted self-generated flux 
4> s (x) has a typical amplitude of variations .which is much 
less than the one observed experimentally£j 



the total phase difference ip(x). To solve this contradic- 
tion we write as 



simp 



[G{x) - G a (x)} . 



(31) 



The function G a (x) is a smoothing average of G(x) over 
an interval with a certain length C ai where I < C a *C C 
The procedure of filtering out the smooth part of G{x) 
is especially evident if we use Fourier series for g{x) and 
G{x). Introducing Fourier transform for g(x) as 



9k e 



(32) 



we find Fourier series for the function G(x) in the form 

00 

G (*) = -E§ efe ' ( 33 ) 



IV. NON-PERIODIC ALTERNATING CRITICAL 
CURRENT DENSITY 

The above analytical approach to the problem of a self- 
generated flux in a non-uniform Josephson junction is 
based on an assumption that the critical current den- 
sity j c (x) is an alternating periodic function. This model 
allows for analytical calculation and provides a reason- 
able preliminary insight into the properties of an ideal- 
ized Josephson junction with an alternating j c {x). At 
the same time this simple model fails for a quantitative 
description of any real system with a certain randomness 
of the spatial distribution of Josephson critical current 
density j c (x). In this section we generalize the above ap- 
proach assuming that the alternating critical current den- 
sity j c (x) is almost periodic, i.e., we assume that there 
is a typical length of interchange of sign of j c {x) which is 
distributed randomly near some mean value 

In the case of an almost periodic j c (x) we can not ap- 
ply the two-scale perturbation theory in the same way as 
we did it in the previous section. Indeed, an arbitrary 
solution of Eq. (p) takes the form 



' ^2 C( x )' 



where 



G(x) = / dx' dx" g{x") 



(29) 



(30) 



and the integration constants a and a 1 are defined by 
the boundary conditions. The random function G(x) in- 
creases with an increase of the integration interval. In 
general, the value of \G(x) | can become arbitrarily large if 
the length of the tunnel junction C becomes large enough. 
This is in a contradiction with our main assumption that 
the phase £(x) is a small and fast varying component of 



where the sums in Eqs. ( |32|) and ( p3[) are over the wave 
vectors k — 2nn/C and n is an integer. The smooth part 
of the function G{x) can be obtained by extracting the 
fast Fourier harmonics, i.e., by extracting from the sum 
( |33| ) terms with wave vectors \k\ > k a ~ 2Ti/C a . As a 
result we find for G a (x) and £(x) the series 



G a (x) 



ka 

-E 


9k ikx 


~k a 




simp 
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(34) 
(35) 



The small and fast alternating part £(x) of the phase 
difference tp{x) is thus defined by Eq. (|3l|). This equa- 
tion is a straightforward generalization of the two-scale 
perturbation theory approach to a real case of an almost 
periodic g{x). 

Next, we use Eq. (|l|) to derive an equation for the 
smoo th p art tp(x) of the phase <p(x). First, combining 
Eqs. (|3l"|), (^), and ( |30| ) we arrive to the relation 

Ajip" - [1 + Gl{x)] simp - g{x)£,{x) sin V = 0. (36) 

Second, we average Eq. ([36]) over an interval with a cer- 
tain length L assuming that I -C L <C C a . This averaging 
results in an equation describing the phase ip(x) 

A 2 jtp" - [l + Gl{x)\ sin ip + j(x) sin V> cos V> = 0, (37) 

where 7(2;) is defined by Eq. (|l5|). It is worth mentioning 
that equation ( |37| ) differs from the analogues equation 
@ by an additional term G'^(x) in the coefficient before 
sin-0 and by the fact that the parameter 7 = j(x) is a 
function of the coordinate x along the junction. 

The coefficient 1 + G'^(x) is defined by magnetic field 
B c (x) which would be generated by a current with the 
density j c (x). Indeed, using Eqs. ([!]), (J30|) , and Maxwell's 
equation dB c /dx — Air j c /c we obtain for G(x) 



4 



G(x) 



dx' [B c (x') - B a (x')] = 



^[<f, c ( x )-ct, a (x)}, 



(38) 



where the magnetic field B a (x) would be generated by 
a constant current density (j c ), i.e., dB a /dx = 4ir(j c )/c, 
and 4> c {x) and 4> a {x) are the fluxes of the fields B c (x) and 
B a {x). It follows now from Eq. (|38|) that 



l + G'i(x) = H^i(^(x))' 

00 



(39) 



As shown above the value of the parameter 7 is de- 
termining the existence or absence of the self-generated 
flux. Using Eq. (||) and the definition of 7(2) given by 
Eq. (plf) we obtain for 7(2) an expression 



sin " " " 



A2/ (C' 2 )>0 (40) 



sin ip 

demonstrating that the condition 7(2;) > holds also in 
the case of an almost periodic critical current density. 



A. Self-generated flux in a stationary state 



The energy of a Josephson junction £ ip given by 
Eq. ( |22| ) can be written in terms of the smooth part of the 
phase difference ip(x). In the case of an almost periodic 
critical current density j c {x) this equation reads 



f — 



4e 

AjV' 2 - 2cos-0[l + G"(a 



7(2;) sin' 



yj] (43) 



and if 7 r (^) < 1, then the energy £ v has a minimum for 
if) — and a maximum for ip — tt. 

In the regions where the function j r {x) > 1 equation 
( [4l| ) has four solutions ip — —ipj(x), 0, ip~ ( (x), ir, with 



ipJx) 



lr{x) 



1 + G»{x) 



(44) 



The energy £ v has a minimum for ip = ± ip 7 (x) and a 
maximum for -0 = 0, tt. The self-generated flux is thus 
non-zero in the regions with j r > 1. This flux has a fast 
and a smoothly varying parts defined by £(x) and 7 (x). 



Let us now consider stationary solutions for the smooth 
part ip{x) of the phase difference <p(x) qualitatively. As- 
sume first that there are sufficiently large intervals with 
lengths Li 3> Aj, where the function ift(x) is constant or 
varies with a typical space-scale of order Li ^> Aj. In 
this case equation (|37|) reduces to 



[l + G" a {x) - 7(x)cos0] sin-0 = 0. 



(41) 




71/2 




FIG. 2. The dependence of 1/7 on 7. 



FIG. 1. Scheme of the phase difference distribution ip(x) 
as an alternation of two types of solutions corresponding to 
different values of 7,-. 

This equation is similar to Eq. (]24|) which we derived 
for the case of a periodic critical current density j c (x) 
and has different solutions depending on the value of the 
parameter 

In the regions with r ) r {x) < 1 equation ( f4l| ) has two so- 
lutions ip = and ip = 7T and therefore as it follows from 
Eq. ( p9| ) there is no self-generated flux in these regions. 



The randomness of the function g(x) causes variation 
of 7r(^) along the junction. As a result of this varia- 
tion intervals with 7r(x) > 1 are interlaced with inter- 
vals with 7r(x) < 1- As it was mentioned above, in 
the case of -y r > 1 the energy of the Josephson junc- 
tion £ v has a minimum for t/j = ± ip 1 [x) and a maximum 
for ip — 0, 7T. When the value of 77.(2:) changes from 
j r {x) > 1 to j r (x) < 1 the energy £ v still has a maxi- 
mum if ip = 7T, but a state with ip = becomes a state 
with a minimum energy. 

The above results provide a qualitative description 
of experimentally observable flux distribution along a 
Josephson junction with an almost periodic alternating 
critical current density. This flux distribution spatially 
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averaged by the measurement tools is denned by the func- 
tion ip(x) (see Fig. 0). Inside the intervals with 7, (2;) > 1 
the phase ipix) tends to one of the solutions ± tp 7 (x). The 
profile of the function ip(x) correlates with the profile of 
tp 1 {x), though does not coincide with it exactly because 
the solution ip(x) — ?/> 7 (x) was obtained under assump- 
tion i/j" = 0, which does not hold exactly for the intervals 
7 r (x) > 1. The smooth part of the phase difference inside 
the intervals with "f r (x) < 1 is tp = which is consistent 
with the assumption ip" — 0. 

The value of tp^ increases quite fast with an increase of 
the parameter 7 (see Fig. ||). In particular, for 7 = 2 the 
value of is already about 0.75 of its maximum value 
7r/2. This means that for most of the experimentally ob- 
servable peaks of the self-generated flux the values of ip 
will be close to 7r/2 which corresponds to a magnetic flux 
cj>0 /4. In some places of the junction the phase ip changes 
from —ip-y to -07- The flux localized in this area of the 
junction will be close to <j>a/2. 



V. NUMERICAL SIMULATIONS 

A. The finite difference scheme 

To study the self-generation of magnetic flux in a tun- 
nel junction with an alternating critical current den- 
sity numerically we introduced time dependence into the 
main equation (H) 



+ aip- AjV" + [1 + g(x)] sin if = 0, 



(45) 



where a ~ 1 is a decay constant. This approach allows 
to study both dynamics and statics of the system. 

The term atp introduces dissipation. As a result of 
this dissipation the relaxation of the system ends up in 
one of the stable stationary states described by a cer- 
tain solution of the static equation (|4|). Moreover, for 
a given distribution of the critical current density j c (x) 
we obtained different static solutions when we start the 
numeric simulation from different initial states. We com- 
pare and classify these solutions based on the features of 
the function j c (x). Indeed, this function essentially de- 
scribes the pinning properties of the junction. Therefore 
a variety of initial conditions can converge to a similar 
flux pinning pattern. 

To solve Eq, ( [45| ) numerically we use the finite differ- 
ence scheme.BIl We adopted this method to our case and 
checked stability and convergency of the obtained solu- 
tions. As a result we arrived to the following scheme 



Yr, 



-Wr, 



Yr, 



-1 - w 



h 2 



(46) 
(47) 
(48) 
(49) 



where /" = f(x m ,t n ), r and h are steps along t and x 
correspondingly. Next, we choose units providing A,/ = 1 
and set h=r. As a result we arrive to the following finite 
difference scheme 



Yn 



- 1 = _(i-ar)^- 1 + (2-arr)^ 



! (1 - g m )smip 



in ■ 



(50) 



B. Stationary solutions 

Initially a certain random function g{x) is generated 
for an interval with a length L with a given values of I 
and 81 (a typical length-scale of the function g(x) and 
its dispersion), g and Sg (amplitude of the function g{x) 
and its dispersion). This allows to calculate the function 
j r (x) for the whole interval. An initial state tpo(x) is 
prepared as a random or some specific function. Finally 
the dynamical rules (j5(]) are applied to the initial state 
itcrativcly until a stationary state is established. 

In Fig. U we show one of the stationary solutions ob- 
tained by a numerical simulation and the function 7r(^_) 
calculated for this solution. It is clearly seen from Fig. y 
that ip(x) arises at the places where 7 r (z) exceeds 1. 
Heights of the peaks are less than 7r/2, and thus the cor- 
responding magnetic flux amplitudes are less than </>o/4. 
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FIG. 3. Profiles of y r {x) and the corresponding ip(x). 

In general a different initial state of the same sample, 
i.e., for the same function 7 r (x), generates a different 
stationary state. Our numerical simulations show that 
these different states differ only by sign of some peaks of 
ip(x), but the shapes and locations stay unchanged. 

WeJaave compared our results with the experimental 
data£j The typical amplitude of the flux variations mea- 
sured by a SQUID pickup loop with a size of several Aj 
is about 0.25 of </>o with rare narrow picks with an am- 
plitude about 0.5 4>o which is in a good agreement with 
our calculations. 
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VI. SUMMARY 



We treated a Josephson junction with an alternat- 
ing critical current density j c {x) as a model for con- 
sidering electromagnetic properties of grain boundaries 
in YBa2Cu307_ x superconducting films. The study is 
mainly focused on a specific case of an almost periodi- 
cally alternating function j c (x). We demonstrated both 
analytically and numerically that under certain condi- 
tions a self-generated flux pattern arise for this type of 
spatial distribution of the critical current density j c (x). 
The obtained flux pattern with two types of interlacing 
flux domains is similar to the one which was recently ob- 
served experimentally in YBa2Cu307_ a; superconducting 
films in the absence of an external magnetic field. The 
typical amplitude for the magnetic flux peaks is of or- 
der </>o/4 and 4>q/2 and the typical distance between the 
peaks depends on the spatial distribution of j c (x). 
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